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Abstract. Energies of a certain class of fuUerene molecules (elongated, contracted, 
and regular icosahedral fuUerenes) are numerically calculated using a microscopic 
description of carbon-carbon bonding. It is shown how these results can be interpreted 
and comprehended using the theory of elasticity that describes bending of a graphene 
plane. Detailed studies of a wide variety of structures constructed by application of 
the same general principle are performed, and analytical expressions for energies of 
such structures are derived. Comparison of numerical results with the predictions 
of a simple implementation of elasticity theory confirms the usefulness of the latter 
approach. 
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1. Introduction 

Recent years have been marked by an accented interest of scientists in molecules and 
materials made exclusively of carbon atoms (C). This interest has been initiated by 
discoveries of fuUerenes |T| and carbon nanotubes |2], and is fueled by many different 
nanoscopic shapes and structures made of carbon atoms that are continually appearing 
either as a result of experimental studies |31I11IS] or imaginative theoretical constructions 
ini m El in] • structures of interest to this work are those in which carbon atoms are sp^ 
hybridized, as in the graphene layers in graphite. Fullerenes and carbon nanotubes fit 
in this category, but there are many other different shapes and objects that can be 
constructed. The multitude of possible shapes is a result of pronounced anisotropy 
of carbon-carbon interactions and the relatively low energy needed to transform a 
hexagonal carbon ring into a pentagonal one. Inclusion of heptagonal rings allows for 
creation of additional, negative-Gaussian-curvature shapes jH] that are not a subject of 
this article. Carbon shapes are of interest since the same elementary constituents (sp^ 
hybridized carbon atoms) can be "assembled" in a variety of stable shapes that will 
not spontaneously decay or transform in some other shapes once created. The resulting 
shapes and structures may differ in their functionality, whatever their purpose may be. 
The assembly of simple elementary pieces into a possibly "engineered" or "designed" 
shape is a dream of nanotechnologists, and the carbon structures are thus the ideal 
benchmark system. Nanoscale shapes of biological interest are often self-assembled 
- a prime example is the selfassembly of viral coatings from individual proteins that 
make it ^Oj- Even more intriguingly, the symmetries and shapes of icosahedral viruses 
are directly related to symmetries of "gigantic" fullerenes [TH [T21 CB] • Both types of 
structures are characterized by pentagonal "disclinations" Jl] (pentagonal carbon rings 
in fullerenes fl] vs. pentameric protein aggregates in viruses ^H]) in crystalline sheet, 
and both can be constructed from the triangulation of an icosahedron. This may be a 
coincidence, but may also point to more general and fundamental principles that may 
still be discovered in the study of carbon structures and shapes. Viruses, as carbon 
shapes, also show the property of multiformity. The same viral coating proteins may 
assemble in icosahedral shells of varying symmetry, but may also make spherocylindrical, 
conical particles or even protein sheets, depending on the conditions ^3 . Similar 
features are found in sp^ carbon shapes. Although various shapes can be, more or 
less easily imagined, this does not mean that they should have large binding energies, 
or that they should be bound at all. Therefore, it seems pragmatic to have a rule 
of thumb estimate for their energy, based solely on their imagined shape. Theory of 
elasticity provides an excellent framework for such estimates and has been in the past 
successfully applied for this purpose jI21 CHI HEl HI] • The aim of this article is to show 
how the energies of the various shapes that can be imagined may be estimated from 
simple relationships that are quite accurate when applied properly. In certain aspects. 
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this article is a continuation of efforts started in Refs. ^Hl and ^7j, but differs from 
them in that it studies a certain class of convex carbon shapes that include both capped 
carbon nanotubes and fuUerenes as special limiting cases. It also applies the theory of 
elasticity in more difficult circumstances in which the symmetry of shapes is reduced. 

In Sec. El I shall present the geometric construction of a certain class of shapes 
that is of interest to this article. In Sec. Ell shall describe the procedure that is used 
to find minimal energy shapes based on the geometric construction discussed in Sec. 
121 The procedure is based on the implementation of conjugate gradient technique |2I| 
in combination with the latest Brenner's potential for the description of carbon-carbon 
bonding [221 , which provides an excellent opportunity to further examine the predictions 
of this relatively recently proposed potential, although the main message of this article 
is totally independent on the form of the potential used. Section jH contains numerical 
results of this article, together with elements of elasticity theory that are applied in the 
analysis of these results. It will be demonstrated how a simple characterization of the 
graphene elasticity together with the knowledge of the energy for creation of pentagonal 
disclination enables a reliable estimate of energies of various convex shapes. Section |21 
summarizes and concludes the article. 

2. Construction of shapes of interest 

The shapes of interest to this work can be constructed as depicted in Fig. ^ The 
procedure is to cut out the outlined shape from the graphene plane and fold it 
into a polyhedron. The shapes that shall be considered can be loosely termed as 
elongated, regular or contracted icosahedra. Note that the construction results in twelve 
pentagonal carbon rings situated in vertices of (contracted or elongated) icosahedron. 
The equilibrated, minimal energy shapes obtained from the pieces of cut-out graphene 
plane are also displayed in Fig. The numerical procedure that was used to obtain 
these shapes will be discussed in detail in Sec. El At this point, not only how the 
elongated shapes resemble carbon nanotubes, while the highly contracted shapes look 
almost like two cones glued together at their bases - there are two pentagonal carbon 
rings at the "poles" of the shape (upper and lower points in shape depicted in panel 
(b) of Fig. P), and the remaining ten pentagonal carbon rings are arranged around an 
"equator" of the shape, in the polygon vertices of two parallel large pentagons rotated 
with respect to each other by an angle of 7r/5 (the vertices of one of those pentagons 
are denoted by points Pi,...,P5 in Fig. (H). The shapes can be uniquely characterized by 
the lengths of two vectors, Ai and B, since |Ai| = IA2I (see Fig. P). The lengths are 
integer multiples of the distance between the centers of neighbouring carbon hexagonal 
rings, i.e. 
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Figure 1. Geometric construction of tlie "armchair" shapes of interest. The shape 
is cut out from the graphene plane and folded in an elongated (shape in panel (a)) 
or contracted (shape in panel (b)) or regular (not shown) icosahedron. Note that the 
points Pi, P5 (panel (b)) are situated in polygon vertices of a (large) pentagon once 
the shape is folded. 

where a is the nearest neighbor C-C distance in graphite. I shall adopt the convention 
that negative values of n correspond to contracted shapes, so that the shapes in Fig. 
^ can be described by (m = 2, n = 2) (elongated icosahedron) and (m = 2,n = —2) 
(contracted icosahedron). The (regular) icosahedral fuUerenes are obtained for n = 0. 
The condition that has to be fulfilled by m and n integers is n > — m. For n ^ m, 
the construction results in capped armchair (5m, 5m) single-wall carbon nanotubes and 
that is the reason for giving the attribute "armchair" to the chosen subset of shapes. 
The total number of carbon atoms (A^) in the constructed shapes is 

N = eOm^ + 20mn. (2) 

Similar geometrical construction can be used to generate to "zig-zag" or chiral (helical) 
shapes jlH] . However, as the main intention of this work is in application of a continuum 
elasticity theory, the precise structure and atomic symmetry of the shape is of no 
importance, and only its shape matters in this respect. It should be noted here that the 
construction results in single-shell or single-wall caihoia shapes/molecules. All the effects 
that are specific to multishell structures, in particular van der Waals interaction between 
the different shells are not treated in this work. For the application of the continuum 
elasticity theory in presence of van der Waals interactions and/or surrounding elastic 
medium see e.g. Refs. [T?H 1^. 
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3. Finding the minimum-energy shapes 

The mathematically constructed folded polyhedra will not represent the minimum 
energy shapes. Carbon atoms in a minimum-energy structure should be relaxed so 
that the nearest-neighbour C-C distances are close to their optimum value energywise. 
To account for these effects, one should have a reliable description of the energetics of 
carbon bonding. In this work, the C-C interactions are modeled using a relatively recent 
second-generation reactive empirical bond order potential by Brenner et al [22] ■ In this 
model, the potential energy (Eh) of a carbon structure is given by 



where r^j is the distance between nearest-neighbour atoms i and j, and are 
pair-additive repulsive and attractive attractions, respectively, and bij is the bond order 
between the atoms i and j that approximately accounts for the anisotropy of C-C 
bonding and for the many-body character of C-C interactions by including effectively 
the three-body contributions to Et,. The detailed description of the potential, together 
with all relevant parameters can be found in Ref. |22]- Similar potential models have 
been successfully employed in previous research of carbon structures (see e.g. Refs. 



The geometrically constructed, folded shape is used as the initial guess of the 
structure. The local minimum of Ef, in the configurational space spanned by coordinates 
of all carbon atoms is found using the conjugate gradient technique described in Ref. 
fn\ . The conjugate gradient technique optimizes all the coordinates at once, proceeding 
through a sequence of steps along the "noninterfering" or "conjugate" direction in the 
configurational space, so that a minimization along a particular direction does not spoil 
the effect of minimization in other conjugate directions. Similar procedures have been 
used for the same purposes in Refs. [121 ITH| ITl]. 

4. Energies of shapes and considerations based on elasticity theory 

4.I. Infinitely long carbon cylinders (non-capped carbon nanotubes) 

Conceptually the simplest shape that can be constructed from the graphene plane is the 
open-ended (non-capped) carbon nanotube. It does not result from the construction 
described in Sec. 121 but shall be examined first to enable a simple insight in the 
application of elasticity theory to the shapes of interest. I shall consider infinitely long 
single-walled armchair carbon nanotubes. The energetics of such shapes can be studied 
by applying the (one-dimensional) periodic boundary conditions to the carbon-nanotube 
unit cell (ring or several rings of atoms), so that the effective number of independent 
coordinates that have to be optimized with respect to total energy of the nanotube is 
relatively small. Figure |21 displays the energy of infinite, armchair- type carbon cylinders 
as a function of the squared cylinder curvature given by k = 1/ (Rcyi)- The mean radius 




(3) 
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of the cylinder, (Rcyi) was calculated from the relaxed atomic coordinates as 



{ncyi} - — , (4j 

where N is the total number of atoms in the unit cell, rj = {xi,yi,Zi) is the position 
vector of the i-th carbon atom, and the cylinder axis coincides with the z-axis. The 
cylinder radius calculated from Eq. (0} is always very close to the analytical prediction 
for tube radius, R = a/{2 sin[7r/(3m)]}. For example, for (10,10) carbon nanotubes, 
Eq. (jH) and the numerically calculated values of yield nanotube radius of 6.802 
A, while the analytical prediction is 6.794 A. The agreement becomes slightly worse 
with the decrease in the nanotube radius. The energies displayed in Fig. |21are in fact 
excess energies per unit area of the structure (denoted by Ac; Ac = 3a^-\/3/4) and were 
calculated as 

j rpper atom Tpper atoin\ ^ fr\ 

\ cylinder graphene) r^^^^j"^": \"/ 

where E^yl^^^^ and E^glaphlne energies per atom in the cylinder and in the 

graphene plane, respectively. The calculations were carried out for two versions of 
the carbon-carbon potential which obviously predict different slopes of the dependence 
of the excess energy on squared curvature, but note that both curves are linear up to 
very large curvatures i.e. small radii of carbon nanotubes. Calculations with the early 
Tersoff potential [SB] predict that somewhat larger energies are required to bend a piece 
of graphene plane in a cylinder with respect to those obtained by using the latest Brenner 
potential [22j. For the energy per unit atom in the graphene plane the two potentials 
predict -Eg^Ip^heJI^ = —7.39552 eV (Tersoff's potential) and -Eg^aph*"^ = —7.39494 eV 
(Brenner's potential). More important difference between the two potentials is in the 
equilibrium C-C bond lengths they predict [a = 1.4605 A and a = 1.4204 A for the 
Tersoff and Brenner potential, respectively). The slopes of the lines shown in Fig. |21 
(co) were obtained from fitting the calculated results pertaining to (5, 5), (6, 6),.. .,(16, 16) 
carbon nanotubes, and for the Brenner and the Tersoff potential they are 0.863 ± 0.003 
eV and 1.093 ± 0.006 eV, respectively (Tersoff obtained Cq = 1.02 eV in Ref. HH]). 

The demonstrated linear dependence of the excess energy per unit area on squared 
(mean) curvature can be easily understood from the examination of the geometry of bent 
network of graphene bonds. Consider a piece of graphene plane bent on a large cylinder 
of radius a. The most important change in the energetics of carbon network is due 
to a change in angles between the bonds, while the bond lengths remain the same as in 
graphene. Only the bond order term, bij depends on the angles between bonds (see Eq. 

and Refs. [221 123 )• By purely geometrical considerations it is easy to derive how the 
angles between the bonds change when the graphene plane is rolled onto a cylinder. If 
one assumes the plane to be inextensional, and considers large cylinder radii (a/i? -C 1) 
the change in energy per i-th atom is given by 

(6) 
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Figure 2. Excess energy per unit area as a function of squared cylinder (armchair 
carbon nanotube) curvature for the Tersoff potential 23 (circles) and the Brenner 
potential (squares). The results for (4,4) and (10,10) armchair carbon nanotubes 
in this graph are denoted. 



which can also be written as 
AE _ 1 
'A~ ~ 2 



:V''(a) 



dh 
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(7) 

^i^l 27r/3 

where j is the index of one of the three nearest neighbours of z-th atom, k ^ j is any 
of the remaining two neighbors of i, and the angle between the i — j and i — k bonds 
is 9ijk I Explicitly written factor of 1/2 is a consequence of the fact that the energy in 
the Tersoff-Brenner potentials is situated in C-C bonds, and that each bond is shared 
by two atoms. Equation ((Zj) is most easily derived by considering armchair or zig-zag 
carbon nanotubes, since in these cases one of the bonds is either parallel (zig-zag) or 
perpendicular (armchair) to the cylinder axis, but for a/R -C 1, the relation holds also 

X The bond order function bij includes summation over k, the remaining two neighbors of i, k ^ j (see 
Refs. |22l I23p . When evaluating the derivative in Eq. lO, all the indices {i,j,k) are to be considered 
as fixed, but the summation over k in bij of course remains. 
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for helical nanotubes, i.e. for small curvatures the bending energies do not depend on 
the symmetry of the carbon nanotube. Note again that the analytical expression (j7|) 
accounts for the excess energy only in the limit of small curvatures k, i.e. for radii of 
curvature that are large with respect to the graphene lattice constant. 

Equation ((Zj) shows that the excess energy per unit area is proportional to the 
squared curvature, in accordance with the results produced by calculations displayed in 
Fig. 121 (note, however, that only a curvature along one direction is accounted for by Eq. 
dZj) which is sufficient for cylindrical surface). It furthermore predicts that the slope of 
this dependence is given by the properties of the interatomic potential (intriguingly, the 
factor of proportionality does not depend on the form of the repulsive part of potential 
for the class of bond order potentials given in Refs. fl^ 1^). The slope Cq can be 
easiliy read out from Eq. ((Zj), and is calculated to be cq = 0.8301 eV and cq = 1.0152 
eV for Brenner's and Terssof 's potentials respectively. This is fairly close to the slopes 
obtained numerically. The small difference can be attributed to the fact that Eq. ((7j) is 
valid only for small curvatures, (or large radii, a/R -C 1), i.e. it includes only the lowest 
order contribution of curvature (quadratic). For tubes of small radii, higher order terms 
are expected to be of importance. In fact, if one restricts the linear fit in Fig. El to 
smaller curvatures (from (12,12) to (16,16) nanotubes, a/R < 0.175) one finds that the 
slopes are 0.839 eV and 1.036 eV, which is in significantly better agreement with the 
prediction of Eq. Note that for (5,5) nanotubes, a/R = 0.414, which means that 
the tubes of such small radii are already in the region where condition a/i? -C 1 does 
not hold, and the predictions of Eq. are not expected to be very accurate. 

The key conclusion of this subsection is that the excess energy per carbon atom 
(AE") of an infinitely long carbon cylinder can be reliably estimated from 

where Cq is the elastic constant of the graphene plane related to the energetics of its 
bending (the bending rigidity), which, for infinitensimal curvatures, can be calculated 
from the knowledge of the interatomic interactions in graphene as in Eq. ([7|). This 
equation holds irrespectively of the details of the interatomic interactions, which only 
change the value of Cq but not the functional dependence displayed in Eq. (jSl), which was 
illustrated by the examination of energetics predicted by two different models of C-C 
interactions [221123 • Equation (jH]) can also be understood without reference to atomic 
structure of the carbon nanotube by interpreting its left hand side as the bending energy 
per unit area of the material. 



4-2. Icosahedral fullerenes 

Icosahedral fullerenes are (m, n = 0) subset of shapes considered in Sec. |21 They are 
more complex than the infinite carbon nanotubes from the standpoint of elasticity theory 
since they contain twelve pentagonal disclinations situated at vertices of an icosahedron. 
The excess energy of icosahedral fuUerene relative to a piece of infinite graphene sheet 
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with the same number of atoms is thus a result of both bending (and possibly stretching) 
of the sheet and the energy required to create disclinations. In order to calculate the 
bending contribution to the energy, one should have some information on the geometry of 
the minimal energy shape of the icosahedral fullerene. The pentagonal disclination can 
be constructed in a graphene plane, starting from the center of a hexagon, and cutting 
the plane in two directions that make an angle of 60 degrees - these cuts can be easily 
seen in the geometrical constructions in Fig. ^ The smaller, cut-out piece of graphene is 
discarded and the remaining (larger) piece of graphene is then folded to rejoin the edges 
of the cut. The shape of minimal energy that the remaining piece of graphene adopts 
is a cone, since it costs less energy to bend the C-C bonds than to stretch them (it is 
convenient to imagine a circular piece of graphene plane with a hexagon in its center). 
The shape will in fact be a truncated cone (or a conical frustum) , since the pentagonal 
ring at the top of the cone will be flat and parallel to the larger base of the cone. 
The precise description of the shape depends on the values of two-dimensional Young's 
modulus and bending rigidity. For materials like graphene (whose stretching requires 
much more energy than bending), the conical frustum is an excellent approximation of 
the minimal-energy shape. From purely geometrical considerations, the half angle of the 
cone is easily shown to be sin~^(5/6). Shapes of crystalline, hexagonally coordinated 
membranes with pentagonal disclinations have been studied in Ref. and many of 
the features discussed there apply also to the case of graphene. The icosahedral fullerene 
is a bit more complicated from a single cone since it contains twelve pentagonal carbon 
rings. It seems plausible that the fullerene shape should be well described by union of 
twelve truncated cones. The radius of the larger base of a particular cone is restricted 
by the presence of other cones. To illustrate these geometrical considerations, in Fig. El 
the shape of (m = 4, n = 0) icosahedral fullerene is displayed (left column of images), 
which is in a minimum of total energy, obtained numerically as discussed in Sec. 0] The 
images in the right column of Fig. IHl display a union of twelve cones with half angles 
of sin~^(5/6) whose larger bases just touch. Note that this geometrical construction is 
an excellent approximation to the calculated minimal energy shape. Its small drawback 
is that its total area (without inclusion of the areas of the larger bases of the cones) 
is smaller from the total area of the fullerene shape which is manifested is Fig. El by 
the appearance of "holes". Note also that the regions where these holes appear are 
almost flat in the numerically calculated minimal energy shape of the fullerene, and 
their contribution to the bending energy should thus be small. 

Having such a good representation of the shape, it is possible to directly proceed 
to the calculation of its bending energy. In the continuum elastic theory of membranes 

the bending energy (Hb) is represented as 



where Cq is the bending rigidity discussed in the previous subsection, cg is the Gaussian 
rigidity, dS is the infinitensimal element of the shape surface, and k and 7 are the mean 
and Gaussian curvature of the surface, respectively (if Ri and R2 are the principal radii 




(9) 
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Figure 3. Panel (a): Equilibrated, minimal energy {m = 4:,n — 0) shape (left) and 
its approximation by twelve conical frusta whose smaller bases centers are situated at 
the vertices of an icosahedron (right). Panel (b): Same as in panel (a), only viewed 
from a different point. Panel (c): Elementary building block of the shapes shown in 
right columns of panels (a) and (b) - conical frustum whose half-angle is sin" ^(5/6). 
The quantities relevant to Eq. Hll|l are denoted. 



of curvature, n = R^^ + R2^, and 7 = {RiR2)~^)- Note that for a cylindrical surface, 
7 = and k = l/Rcyi {Ri = Rcyi, R2 = cc) in every point of the surface, and the 
total bending energy integrates to Hh/A = where A is the total area of the 

cylinder (without the bases). This is in agreement with Eq. (jHJ. For conical surfaces, 
the gaussian curvatures are also zero everywhere, so the total bending energy of the 
approximated fuUerene shape will depend only on cq. Integrating Eq. over the 
surface of a cone whose half-angle is /3, one finds that the bending energy of a cone {H^) 
is 

cos^ P 

= . n In {rmax/rmin) , (10) 

sm p 

which for the cones with P = sin~^(5/6) yields [TE] 
llvr 

-f^b = "20"*^° ^n i^max/f^min) ; (H) 

where the meaning of and Vmax is illustrated in panel (c) of Fig. El - these are 
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the distances of the upper and lower bases of the conical frustum from the cone apex 
measured along the cone face. The value of r^m can be easily evaluated due to the 
fact that the smaller base of the frustum is a pentagon of carbon atoms. There is 
some ambiguity, however, since the pentagon is to be approximated by a circle. If the 
circumradius of the pentagon is identified with the radius of the smaller base of the 



cone, one obtains that r^m = Say 50 + 10-\/5/25 ~ 1.021a, which is the upper bound 
for Tmin- The identification of the radius of the smaller base of the cone with the inradius 



of the pentagon yields rmin = Say 25 + 10-\/5/25 ~ 0.826a, which is the lower bound 
for Vmin- In any case, Tmin = fa, where / is a numerical factor between 0.826 and 
1.021. The value of r^ax is half of the shortest distance between the two apices of the 
neighboring cones measured along the cone faces. From Figs. Q and El one finds that 
^max = 3(m — 2)a/2 + r^m, and by using Eq. Q one obtains that 




(12) 

60 / ^ ^ 

It is obviously quite convenient to choose / = 1, closer to the upper bound for r^m, so 
that the total excess energy of the icosahedral fuUerene can be written as 



AE{N) = 12 



llvr , /A^\ llvr , , . 

M H Cn In — H Cn In - 

. 60 ° V60/ 30 ° * " 



+ Eholes- (13) 



The multiplicative factor of 12 in the above equation is due to sum of the energies of 
twelve cones that make the shape surface, and A5 is the core energy of the pentagonal 
disclination which contains the local effects associated with the atomic structure of 
carbon pentagon and its immediate neighborhood. The contribution of the part of the 
fuUerene surface that was not covered by the union of cones to the total energy was 
denoted by Ehoies- Note that this quantity should represent a quite small correction to 
the total energy (see Fig. Oj) and is not expected to depend strongly on A^, since the 
area of the "holes " is proportional to A^, while their mean curvature is approximately 
proportional to (inverse radius of the fuUerene, see Eq. Q). Although one could 

easily invent a nice geometrical model to estimate this contribution, it is possible to 
proceed further by noting that Eq. (fT3|) can also be written as 



AE{N) = AE{Cm) + 



llvr , fN 
Co In 



(14) 



5 V60. 

where A.E{Cqq) is the total excess energy of Ceo which the calculation with the 
Brenner potential predicts /S.E{Cqq) = 26.744 eV. Figure El displays the results of 
numerical minimization of total energy of icosahedral fuUerenes together with the 
prediction of Eq. ()14|) with cq = 0.863 eV taken from the analysis of infinite carbon 
cylinders. It seems more sensible to choose this value of co rather than the one predicted 
by Eq. ^ which is appropriate for small curvatures, since close to the smaller bases 
of the truncated cones, condition a/R <C 1, where R is the radius of curvature is not 
fulfilled. Thus, a usage of a somewhat larger value of Cq obtained from the numerical 
analysis of infinite carbon cylinder should approximately account for the contribution 
of terms like {a/Ry, p > 2 to the bending energy of the shapes studied. Note, however. 



Energies of sp^ carbon shapes with pentagonal disclinations and elasticity theory 12 

that the (infinitensimal) analytical value of Cq (0.83 eV) and the one obtained from the 
study of infinite carbon cylinders (0.863 eV) are quite similar, and one could also chose 
to work with the analytical value of Cq with comparable success (see below). 




ln(A^/60) 

Figure 4. Numerically calculated excess energies of icosahedral fuUerenes (squares) 
and the prediction of Ea. H14|) (line). 

It can be seen in Fig. HI that the agreement between the calculated values and those 
predicted by Eq. (|Tl|l is excellent (note that there are no fit parameters). Very small 
deviations can be mostly attributed to the fact that the regions around the larger bases 
of the cone are constrained by the continuity of the fullerene surface, and deviation from 
the cone shape is thus expected. It is to a certain extent surprising that the predictions 
of Eq. (fT^ are so precise, especially concerning the fact that the excess energy of 
Ceo figures in it as a parameter. Buckminsterfullerene is certainly not expected to be 
reliably described by the elasticity theory, since there are no atoms on the cone faces, 
only C-C bonds. Yet, the elasticity theory results maintain their meaning down to 
fullerenes of very small sizes and small number of atoms. By examining the energies 
of individual atoms in the icosahedral fullerenes it is possible to proceed a bit further 
and to estimate the value of the core energy, A5. For example, in (m = 4, n = 0) 
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fuUerene shown in Fig. |2l the energies of atoms in pentagonal rings are -6.959 eV, while 
atoms situated in hexagonal carbon rings have energies of about -7.38 eV, depending on 
their position (see Eq. (jHl)). The energies of atoms in pentagonal rings do not change 
significantly with increasing the size of the fullerene, which is a consequence of the local 
nature of this energy. All this suggests that A5 = 5 [—6.9587 — (—7.3949)] eV = 2.18 
eV, where -7.3949 eV is the energy per atom in the infinite graphene plane, predicted 
by the Brenner potential. Equation ()14p does not depend (formally) on this quantity, 
but it can be calculated from Eq. (fT!?|) which yields Ag''*** = 1.83 eV {Ehoies was set 
to zero), in a quite nice agreement with the number obtained from the detailed atomic 
description of the icosahedral fullerenes, especially in a view of all the details involved 
in the approximation of the fullerene shape. This also confirms that E^oies is indeed 
negligible part of the total energy. 



AE{m,n) = 12 



. 30 \ 2 J 



+ Eholes, (15) 



4-3. Capped carbon nanotubes, elongated icosahedral fullerenes 

Using the same reasoning outlined in the previous subsection, it is easy to construct 
an approximate shape of the elongated icosahedral fullerene (n 7^ 0). These can be 
approximated by the union of twelve conical frusta and a cylinder whose length is na\/3, 
and radius 15ma/(27r) (a finite-length piece of armchair (5m, 5m) carbon nanotube). 
The excess elastic energy can be thus written as 

27r^Co n 
5^3 m 

which could also be parameterized by the total number of atoms, A^, and n. The first 
term in the expression for energy in Eq. (fTH|) is the sum of the local and bending 
energies associated with the twelve equal pentagonal disclinations (conical frusta), and 
the second one corresponds to the bending energy of the finite length cylinder. The third 
term is again the energy associated with the "holes" in the geometrical construction, as 
discussed in the previous subsection. Note that for n = 0, Eq. (fT3j) reduces to Eq. (fT^ . 
Although the geometrical construction yielding Eq. (|15|) is quite simple, it has a small 
drawback in that the areas of the cones and of the cylinder in between them overlap to a 
small extent. Angles between the surfaces of the cones and the cylinder at the points in 
which the entities nearly touch are small but nonvanishing, which is another drawback 
of the proposed construction. More elaborate constructions are possible, but are not 
needed (Fig. Equation [13 can be approximated by 



AE{m,n) = AE{Q 



60 j 



227r , , , 27i^Co n 
■Co In (m) + 



(16) 



5 5^3 m_ 

which simply states that the excess energy of the shape is equal to the sum of the excess 
energy of the icosahedral fullerene and the bending energy of the cylinder inserted in 
between its two "halves" (see Fig. El). This equation is, however, an approximation, 
since it is based on the assumption that the excess energies associated with the holes 
in the geometrical constructions are the same for the icosahedral fullerene and the 
elongated icosahedral fullerene, which is not exactly the case (compare Figs. El and E}- 
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(a) (b) 

Figure 5. Panel (a): Numerically calculated equilibrium shape of {m = 2,n — 4) 
elongated icosahedral fuUerene. Panel (b): Approximation of shape using geometrical 
construction leading to Eq. H15|l . 

It furthermore contains an assumption that the elastic energies of holes do not depend 
on m and n integers characterizing the structure. Nevertheless, an excellent account 
of the numerical data is obtained which can be seen from Fig. IHl Note again that the 
analytical results presented in Fig. IHl represent the predictions of the elasticity theory 
without fit parameters. 

4.4- Contracted icosahedral fullerenes 

The contracted icosahedral fullerenes are more complicated from the shapes considered 
thus far. Construction of their geometrical approximation is not simple mostly due to 
specific geometrical constraints that the ten pentagonal disclinations situated around 
the equator of the shape are subjected to. When n vanishes, the area around each of 
the disclinations is well represented by a cone with an angle of 2sin~^(5/6) which was 
clearly demonstrated in subsection 14.21 When n decreases (i.e. increases in absolute 
value), the two large pentagons parallel to the shape equator whose vertices coincide 
with five pentagonal disclinations (see Fig. approach each other. This constrains 
the shape of the area around these pentagonal disclinations. As the cone was such an 
useful approximation of the area around pentagonal disclinations for shapes studied in 
subsections 14 . 21 and 14. H| it seems reasonable to try to describe the contracted icosahedral 
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1 2 3 4 5 6 7 8 9 10 11 12 13 

n 




m 

Figure 6. Panel (a): Excess energies of elongated icosahedral fuUerenes (capped 
carbon nanotubes) as a function of elongation number (n) with m as a parameter of 
the curves. Full symbols correspond to numerically calculated excess energies, while 
lines correspond to the predictions of elasticity theory in Eq. (|16|) . Panel (b): Same 
data as in panel (a) with the roles of m and n numbers interchanged. 

fuUerenes also as union of cones, ten (those arranged around the equator of the shape) 
of which belong to one category, and two (the ones at the two poles of the shape) to 
another. The two cones at the poles are easily constructed as in subsections 14 . 21 and 14 . 31 
- these have the angles of 2sin~^(5/6). The area around ten remaining disclinations is 
dominated by two opposing geometrical constraints. As n decreases, the angle subtended 
by a line that connects a particular disclination (point A in panel (b) of Fig. Ej) and 
the one at the closest pole (point B) and a line that connects the disclination with 
the midpoint of a line drawn between the two nearest disclinations in the neigbouring 
pentagonal ring of disclinations (point C) also decreases - this angle is denoted with ai 
in panel (b) of Fig. [7| On the other hand, the angles that are subtended by the lines 
DA and AE or FA and AG (see Fig. Ej) increase when n decreases. It is thus clear that 
the area around the disclination cannot be represented by a single cone. Nevertheless, it 
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seems reasonable to approximate this area by a union of pieces of surfaces of two cones, 
each of them satisfying one of the two opposing geometrical constraints. The angle of 




(a) 




(b) 



Figure 7. Panel (a): Equilibrium shape of (m = 3, n = —2) contracted icosahedral 
fuUcrcne. The characteristic pyramid- like shape is denoted by full lines. Panel (b): 
The pyramid- like shape with denotations discussed in the text. 



one of the cones is about ai, which is given by 

1 



ai{n) 



cos 



-1 



2sin(7r/5) 



cos 



-1 



V^mtan(7r/10) 
3m — 2n 



(17) 



The angle of the second cone (02) that is searched for is larger than ai, and the 
contribution of this half-cone to the total elastic energy of the area surrounding the 
disclination is thus smaller (see Eq. (|Tn|l ). Nevertheless, the second cone is important 
since it contains a part of the total area surrounding the disclination, thus reducing the 
area available to the cone with angle ai. When n = 0, both ai and 02 angles should 
be 2sin~^(5/6). Note that when n decreases, angle 02 changes in a significantly smaller 
range than ai, since it is importantly constrained by the fixed angle between lines DA 
and AH (see Fig. [7|). Thus, I shall neglect the change of 0:2 with n and set it to a 
constant 



a2 = 4sin "'^(5/6) — ai{n = 0) 
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4 sin (5/6) — cos 



-1 



1 



cos~^ [tan(7r/10)] 



sin(7r/5) 



2.04 sin-i(5/6), 



(18) 



which is obviously quite close to 2sin~^(5/6), as it should be, and which was chosen 
so that ^2 + diin = 0) = 4sin~^(5/6), since I am further going to assume that each 
of the two cones takes exactly a half of the total area available to the disclination. 
Together with the choice for a2 in Eq. (fT^ . this assumption shall yield the elastic 
disclination energy in the limit when n = equal to that characteristic of icosahedral 
fullerenes studied in subsection 14 .21 Note that the axes of the two cones do not need to 
be identical - the two pieces of the conical surfaces can be rotated with respect to each 
other at will. It is also irrelevant how the particular cone (half)surface is constructed as 
long as the pieces of it are bounded by the arcs of the two bases subtending the same 
angle from the cone axis and the two lines (along the cone face) that connect the arc ends 
in the upper and lower base. Thus, there is a significant freedom in constructing a shape 
around the disclination, but for all of the thus constructed shapes the buckling energies 
are equal. Note also that for n ^ points D, B, H, and E do not lie in a plane. The 
proposed construction is admittedly rather approximate. For example, one could object 
that the fractions of the total area belonging to each of the two cones should also (at 
least slightly) change with n. Nevertheless, the construction accounts for most important 
geometrical constraints imposed on the area around pentagonal disclinations situated 
around the equator of the shape due to the presence of neighboring disclinations. The 
exact answers to all of the objections regarding the proposed construction that could 
be easily put forth would require a solution of the nonlinear problem in the theory of 
elasticity, which is not the aim of this study - the aim is to obtain a simple insight in the 
energetics of the shapes using the simplest possible application of the elasticity theory, 
and to obtain, preferably analytic, expressions for the excess energies of such shapes. 

Following the reasoning presented in the previous paragraph and using Eq. (fTUI) for 
the calculation of bending energy of a cone with arbitrary half-angle, the total excess 
energy of the proposed geometrical approximation to the contracted icosahedral shapes 
is 



AE{m, n) = I2A5 + 2Epoie{m, n) + IQEequatorijn, n) + Ehoi 



(19) 



where 



E. 



'pole 



llCoTT 

30 



ln(3m — q) 



(20) 





and 



q = 




(22) 



4 
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The last parameter (g) is simply a half of the distance (in units of a) between the two 
nearest neighboring disclinations situated at different pentagonal rings of disclinations 
around the equator of the shape (points A and E in Fig. Ej). 

Numerically calculated excess energies of the shapes and the predictions of Eq. (fT^ 
are shown in Fig. IHland denoted by symbols and lines, respectively . In evaluating Eq. 
finU) . Ehoies was set to zero, A5 = 1.83 eV was taken from the analysis of icosahedral 
fuUerenes in subsection 14.21 and Cq = 0.863 eV, as in the two previous subsections. Note 
again that there are no fitting parameters involved. The agreement of the predictions 




Figure 8. Panel (a): Excess energies of contracted icosahedral fullerenes as a function 
of n with m as a parameter of the curves. Full symbols correspond to numerically 
calculated excess energies, while lines correspond to the predictions of Eq. (|19|l . Panel 
(b): Same data as in panel (a) with the roles of m and n numbers interchanged. 



of the elasticity theory with the numerical results is striking. In spite of all the 
approximations involved in the geometrical construction of the equilibrium shape, Eq. 
(fTTHl both qualitatively and quantitatively accounts for the numerical data. The largest 
disagreement is found for the smallest shapes (m = 2), but the curves correctly account 
for the appearance of minimum in excess energies for some finite n (panel (a) of Fig. |H1), 
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and for the increase in excess energies when \n\ becomes comparable to m. It is worth 
noting that although the number of atoms in n = —m structures is smaller than in 
{m,n = 0) structures (regular icosahedral fullerenes) by 20 mn (see Eq. (j2I)), the total 
excess energies are still significantly larger in such structures, which is a consequence 
of an energetically unfavorable geometrical positions occupied by the ten disclinations 
around the equator of the shape. This consideration can also explain the minimum 
in excess energies obtained for some n when m is fixed - as the elongation number 
n decreases, the number of atoms in the structure decrease, but the contribution to 
excess energy due to ten disclinations around the equator increase. The combined 
consequences of these two effects lead to the appearance of a minimum in excess energy 
as a function of n. Of the three types of shapes studied, the contracted icosahedral 
fullerenes typically have the largest excess energies per atom. For example, the excess 
energy of (m = 10, n = —10) contracted fullerene {AE = 59.47 eV) that contains 4000 
atoms is significantly larger from the excess energy of (m = 9, n = 0) regular icosahedral 
fullerene [AE = 53.16 eV) that contains /ar^'er number of atoms {N = 4860). 

5. Summary and conclusion 

It has been shown that the insight in the geometry of the equilibrium shapes can be 

extremely useful in estimating their (excess) energies. Such an approach does not require 

detailed description of the interatomic interactions characteristic of a shape or molecule 

in question, instead it relies exclusively on the knowledge of elastic parameters of the 

material that the shape (shell) is made of. In the studies presented in this article, the 

parameters that were necessary for reliable estimation of energy were the local energy of 

pentagonal ring of carbon atoms (A5), and the bending rigidity of graphene plane (cq). 

Both of these parameters derive from the quantum-mechanical nature of bonding of 

atoms (see Eq. ((7j)) that make the shape. Surprisingly enough, it has been demonstrated 

that the application of the simplest theory of elasticity yields reliable results even for 

shapes made of a small number of atoms, i.e. in the regime where its applicability is not 

to be expected. For all of the shapes considered, the knowledge of the two-dimensional 

Young's modulus of graphene [21] was not needed, which suggests that in such shapes 

almost all of the elastic energy is of the bending type, while the energies associated 

with stretching are negligible §. For shapes in which the geometrical constraints on the 

shape are strong (as in the case of contracted icosahedral fullerenes), the application of 

elasticity theory may be involved and complicated. Nevertheless, the identification of 

the constraints and their effects on the allowed shapes yields an additional insight in 

the energetics of the shapes and its behaviour within a certain class of shapes studied. 
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